The consequences of nuclear electron capture in core collapse supernovae 
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The most important weak nuclear interaction to the dynamics of stellar core collapse is electron 
capture, primarily on nuclei with masses larger than 60. In prior simulations of core collapse, electron 
capture on these nuclei has been treated in a highly parameterized fashion, if not ignored. With 
realistic treatment of electron capture on heavy nuclei come significant changes in the hydrodynamics 
of core collapse and bounce. We discuss these as well as the ramifications for the post-bounce 
evolution in core collapse supernovae. 

PACS numbers: 97.60.Bw, 25.50.+X, 23.40.-s 



Core collapse supernovae are among the most energetic 
events in the universe, emitting 10 46 J of energy, mostly 
in the form of neutrinos. These explosions mark the end 
of the life of a massive star, the formation of a neutron 
star or black hole and play a preeminent role in the cos- 
mic origin of the elements. With the formation of an iron 
core in a massive star (containing the maximally-bound 
iron and neighboring nuclei) , thermonuclear energy is no 
longer available to slow the inexorable contraction that 
results from a star's self gravity. Once this cold iron core 
grows too massive to be supported by the pressure of de- 
generate electrons, core collapse ensues. In the inner re- 
gion of the core, this collapse is subsonic and homologous, 
while the outer regions collapse supersonically. When 
the inner core exceeds nuclear densities, it stiffens, halt- 
ing the collapse. Collision of the supersonically infalling 
outer core with this stiffened inner core produces the 
bounce shock, which initially drives outward the outer lay- 
ers of the iron core. However, this bounce shock is sapped 
of energy by the escape of neutrinos and nuclear dissoci- 
ation and stalls before it can drive off the envelope of the 
star (see, e.g., The intense neutrino flux, which is 

carrying off the binding energy of the proto-neutron star 
(PNS), heats matter between the neutrinospheres and 
the stalled shock. In the neutrino reheating paradigm, 
this heating reenergizes the shock, which drives off the 
concentric layers of successively lighter elements that lie 
above the iron core, producing the supernova. 

Unfortunately, simulations exploring the neutrino re- 
heating paradigm often fail to produce explosions. The 
failure of recent spherically symmetric multigroup Boltz- 
mann simulations U3- LJ- Lj to produce explosions has re- 
moved incomplete neutrino transport as a potential cause 
of this failure. Models that break the assumption of 
spherical symmetry have achieved some success, cither 



by an increase in the neutrino luminosity due to fluid 
instabilities within the proto-neutron star @ or by en- 
hancement of the efficiency of the neutrino heating by 
large scale convection behind the shock 0, HI @ The 
PNS instabilities are driven by lepton and entropy gradi- 
ents, while convection behind the shock originates from 
gradients in entropy that result from the stalling of the 
shock and grow as the matter is heated from below. How- 
ever, even with such enhancements, explosions are not 
guaranteed 0, 0, 0] . A third potential cause of the 
failure to produce explosions in numerical models is in- 
complete or inaccurate treatment of the wide variety of 
nuclear and weak interaction physics that is important 
to the supernova mechanism. Once the supernova shock 
forms, emission and absorption of electron neutrinos and 
antincutrinos on the dissociation-liberated free nucleons 
are the dominant processes. However, during core col- 
lapse, electron capture on nuclei plays an important and, 
as we will demonstrate, dominant role by significantly 
altering the electron fraction and entropy, thereby deter- 
mining the strength and location of the initial supernova 
shock, as well as the entropy and electron fraction profiles 
throughout the core. As a result, improvements in the 
treatment of electron capture alter the initial conditions 
for the entire postbounce evolution of the supernova. 

Calculation of the rate of electron capture on heavy nu- 
clei in the collapsing core requires two components: the 
appropriate electron capture reaction rates and knowl- 
edge of the nuclear composition. The inclusion of electron 
capture within a multigroup neutrino transport simula- 
tion adds an additional requirement: information about 
the spectra of emitted neutrinos. Unlike stellar evolution 
and supernova nucleosynthesis simulations, wherein the 
nuclear composition is tracked in detail via a reaction 
network Il4| , in simulations of the supernova mecha- 
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nism, the composition in the iron core is calculated by the 
equation of state assuming nuclear statistical equilibrium 
(NSE). Typically, the information on the nuclear compo- 
sition provided by the equation of state is limited to the 
mass fractions of free neutrons and protons, a-particles, 
and the sum of all heavy nuclei, as well as the identity of 
an average heavy nucleus, calculated in the liquid drop 
framewor k Il5l. In most recent supernova simulations 
(see , e.g., I^lsllllt), the treatment introduced by Bruenn 
(17| is used. This prescription treats electron capture on 
heavy nuclei thr oug h a generic O/7/2 — > Of 5/2 Gamow- 
Teller resonance in the average heavy nucleus identi- 
fied by the equation of state. Because this treatment does 
not include additional Gamow- Teller transitions, forbid- 
den transitions, or thermal unblocking (see 0), electron 
capture on heavy nuclei ceases when the neutron number 
of the average nucleus exceeds 40. As a result, electron 
capture on protons dominates the later phases of collapse. 

As a major advance over this simple treatment of nu- 
clear electron capture, we have developed a treatment 
based on recent shell model electron capture rates from 
Langanke & Martinez Pinedo (LMP) |2q| for 45 < A < 
65 and 80 reaction rates from a hybrid shell model- 
RPA calculation [l^. |2~J | for a sample of nuclei with 
A = 66 — 112. For the distribution of emitted neutrinos, 
we use the approximation described by Langanke et al. 
|22^ . To calculate the needed abundances of the heavy 
nuclei, a Saha-like NSE is assumed, including Coulomb 
corrections to the nuclear binding energy I23l|24l , but ne- 
glecting the effects of degenerate nucleons |25j . This NSE 
treatment has been used in prior investigations of elec- 
tron capture in thermonuclear supernovae |2^|. We use 
the combined set of LMP and hybrid model rates to cal- 
culate an average neutrino emissivity per heavy nucleus. 
The full neutrino emissivity is then the product of this 
average and the number density of heavy nuclei calcu- 
lated by the equation of state. With the limited coverage 
of rates for A > 65, this approach provides the most rea- 
sonable estimate of what the total electron capture would 
be if rates for all nuclei were available. This averaging 
approach also makes the rate of electron capture consis- 
tent with the composition returned by the equation of 
state, while minimizing the impact of the limitations of 
our NSE treatment. A more detailed description of our 
method, including tests of some assumptions made, will 
be presented in a forthcoming article 27] . 

Simulations of the collapse, bounce, and post-bounce 
evolution of a 15 M Q model were carried out us- 
ing the fiducial Bruenn prescription for electron cap- 
ture on nuclei and our LMP+hybrid treatment with our 
fully general relativistic, spherically symmetric AGILE- 
BOLTZTRAN code. In these simulations, it is employed 
using the equation of state of Lattimer and Swesty [15j , 
6-point Gaussian quadrature to discretize the neutrino 
angular distributions, and 12 energy groups to discretize 
the neutrino spectra between 3 to 300 MeV. 
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FIG. 1: The electron fraction, entropy, density, and velocity as 
functions of the enclosed mass at the beginning of bounce for a 
15 Mq model. The thin line is a simulation using the Bruenn 
parameterization while the thick line is for a simulation using 
the LMP and hybrid reaction rate sets. 



Our improved treatment of nuclear electron capture 
has two competing effects. In lower density regions, 
where the average nucleus is well below the N — 40 
cutoff of electron capture on heavy nuclei, the Bruenn 
parameterization results in more electron capture than 
the LMP+hybrid treatment. This is similar to the re- 
duction in the amount of electron capture seen in stellar 
evolution models [2^ and thermonuclear supernova j2|| 
models when earlier parameterized rates |29) are replaced 
by shell model calculations. In denser regions, the contin- 
uation of electron capture on heavy nuclei alongside elec- 
tron capture on protons results in more electron capture 
in the LMP+hybrid case. The results of these competing 
effects can be seen in the upper pane of Figure ^ which 
shows the distributions of electron fraction, entropy, den- 
sity, and velocity throughout the core at bounce (maxi- 
mum central density). 

In addition to the marked reduction (~ 10%) in the 
electron fraction in the interior of the PNS, the improved 
treatment of electron capture also results in an ~ 20% 
reduction in the mass of the homologous core, consistent 
with the analysis that the size of the homologous core 
is proportional to the square of the mean trapped lep- 
ton fraction < F; 2 > at core bounce At bounce, 

this change in the homologous core manifests itself as a 
reduction in the mass interior to the formation of the 
shock from 0.57 M Q in the fiducial case to 0.48 M Q in 
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FIG. 2: The neutrino luminosity and root-mean-square en- 
ergy (at 500 km) as a function of time from bounce for a 15 
M© model. The thin lines show this evolution for a simula- 
tion using the Bruenn parameterization while the thick lines 
show this evolution for a simulation using the LMP+hybrid 
rates. The solid lines correspond to electron type neutrinos, 
the dashed lines correspond to electron type antineutrinos. 

the LMP+hybrid case, as is evident in the lower pane 
of Figure ^ In the LMP+hybrid case there is also an 
~ 15% reduction in the central density and an ~ 5% re- 
duction in the central entropy at bounce, as well as an 
~ 15% smaller velocity difference across the shock. 

If the principle effect of our improvement in the treat- 
ment of electron capture on nuclei were to launch a 
weaker shock with more of the iron core overlying it, this 
improvement would make a successful explosion more dif- 
ficult. However, these improvements in nuclear electron 
capture also alter the behavior of the outer layers which 
play an important role in the ultimate fate of the shock. 
The lesser neutronization in the outer layers slows the col- 
lapse of these layers, which further diminishes the growth 
of the electron capture rate by reducing the rate at which 
the density increases. These changes are clearly apparent 
in regions above the shocks in Figure ^ f° r example, re- 
ductions of a factor of 5 in density and 40% in velocity are 
evident in the vicinity of 0.8 Mq. Such changes reduce 
the ram pressure opposing the shock, easing its outward 
progress. In these spherically symmetric models, these 
improvements allow the shock in the LMP+hybrid case 
to reach 168 km, relative to 166 km in the fiducial case, in 
spite of the greater mass overlying the shock when it was 
launched and the greater loss of energy to the neutrino 
burst (see Figure EJ). 

As mentioned earlier, changes in the electron capture 
rates also lead to changes in the core fluid gradients that 
may, in turn, drive fluid instabilities that are potentially 



important to the supernova mechanism. Within the in- 
ner 50 km, the entropy and lepton fraction gradients 
found in the LMP+hybrid model are considerably dif- 
ferent from those found in the fiducial model. Conse- 
quently, the more accurate treatment of electron capture 
may significantly alter the location, extent, and strength 
of proto-neutron star convection, or other potential fluid 
instabilities, in the core. This provides an excellent ex- 
ample of the coupling of convective behavior to radiative 
and nuclear processes and must be further investigated 
in the context of future multidimensional models. 

Figure |21 shows the luminosity and mean energy of 
the emitted electron neutrinos and antineutrinos between 
100 milliseconds before bounce and 100 milliseconds af- 
ter bounce. Clearly evident in the luminosity is a slight 
delay (2 ms) in the prominent "breakout" burst caused 
by the deeper launch of the shock in the LMP+hybrid 
case. Over the first 50 milliseconds after bounce, the 
LMP+hybrid model emits ~ 15% more energy than the 
fiducial model, with a slightly lower luminosity at later 
times. This is largely the result of differences in the mean 
electron neutrino energy, which is as much as 1 MeV 
higher over the first 50 milliseconds in the LMP+hybrid 
case, but lower thereafter. This results from the neutri- 
nospheres in the LMP+hybrid model occurring in deeper, 
hotter layers for the first 50 milliseconds, but cooler lay- 
ers at later times. 

The differences in the neutrino spectrum during col- 
lapse, when electron capture on nuclei dominates, are 
larger than those described after bounce. For low den- 
sities, where capture on nuclei dominates in the Bruenn 
prescription as well, the approximate reaction Q-value 
derived from the free neutron and proton chemical po- 
tentials dramatically underestimates the Q value, result- 
ing in a much lower mean neutrino energy. As captures 
on protons begin to compete with captures on nuclei in 
the Bruenn prescription, the mean neutrino energy grows 
rapidly because of the higher Q value for capture on pro- 
tons. It exceeds that found in our LMP+hybrid model 
by as much as 2 MeV in the 30 milliseconds just before 
bounce. This latter effect was anticipated by Langanke 
et a/.0], who also demonstrated that nuclear electron 
capture should dominate that on protons because the 
much larger number of heavy nuclei more than compen- 
sates for the larger capture rate on free protons, at least 
for the conditions found in models resulting from the 
Bruenn prescription. These self consistent models un- 
equivocally show that this is correct. At the onset of 
collapse, there are roughly 1000 heavy nuclei per proton 
in the inner layers of the core. In our fiducial model this 
ratio declines rapidly, reaching values less than 100 by 
the time the central density is 10 1:L gcm _3 and roughly 
10 by the time the central density is 10 13 gcm~ 3 . In the 
LMP+hybrid model, the reduced electron fraction and 
entropy keep this ratio near 1000 until the central den- 
sity exceeds 10 12 gcm -3 , reaching 50 around a central 
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density of 10 gem". As a result, in the regime where 
Y e experiences the largest changes (10 11-13 g cm -3 ), the 
dominance of heavy nuclei is increased by a factor of 5-30, 
cementing the dominance of nuclear electron capture. 

To implement these simulations, we have made approx- 
imations to all three components of the calculation of 
electron capture in core collapse supernovae. Each of 
these requires further improvement. Removing the aver- 
aging of electron capture rates requires better coverage 
of electron capture on nuclei in the region 65 < A < 
120, by hybrid and other approximate methods [3jJ in 
the near term, but ultimately by shell model calcula- 
tions vetted by experimental determinations of Gamow- 
Teller strength distributions. These reaction rates must 
cover the full thermodynamic range of interest in su- 
pernovae (temperatures of 1-100 GK and densities from 
10 5 — 10 14 gcm -3 ) and must also address the need for 
the emitted neutrino spectra. Detailed tracking of the 
nuclear composition is also necessary, in a form that re- 
tains the consistent transition to nuclear matter afforded 
by current schemes 15] while allowing for accurate cal- 
culation of the rate of electron capture on heavy nuclei 
and, ultimately, for detailed nucleosynthesis. 

We have demonstrated that supernova simulations 
with a modern treatment of electron capture differ sig- 
nificantly from previous models, which employed more 
parameterized treatments. Though this improved model 
still fails to produce an explosion in the spherically sym- 
metric case, the differences are quite striking. The ini- 
tial mass behind the shock when it is launched is re- 
duced by 20%, with significantly (~ 10%) lower central 
densities, entropies, and electron fractions in this region. 
Over the first 50 ms after bounce, the neutrino luminos- 
ity is boosted by ~ 15% with the mean electron neutrino 
energy increased by ~ 1 MeV. In spite of an initially 
weaker and deeper shock and larger neutrino energy loss, 
reduced electron capture in the outer layers slows their 
collapse, allowing the shock to reach a maximum radius 
that is slightly larger. Furthermore the lepton and en- 
tropy gradients in the core differ significantly. Because 
these gradients drive PNS convection and other potential 
instabilities in the core, the location and strength of such 
instabilities may be significantly different than heretofore 
thought. 
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